Na induced correlations in Na x .CoC>2 
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Increasing experimental evidence is building which indicates that signatures of strong correlations 
are present in the Na rich region of Na x Co02 (ie. x « 0.7) and absent in the Na poor region (ie. 
x ~ 0.3). This is unexpected given that NaCo02 is a band insulator and C0O2 has an integer 
filled open shell making it a candidate for strong correlations. We explain these experimental 
observations by presenting a minimal low-energy Hamiltonian for the cobaltates and solving it within 
LDA+DMFT. The Na potential is shown to be a key element in understanding correlations in this 
material. Furthermore, LDA calculations for the realistic Na ordering predict a binary perturbation 
of the Co sites which correlates with the Nai sites (ie. Na sites above/below Co sites). 



The qualitative features of the cobaltate phase dia- 
gram have been experimentally established |l|. Uncor- 
related behavior is observed in the Na poor region while 
correlated behavior is observed in the Na rich region. 
For example, Magnetic susceptibility measurements show 
Pauli-like behavior for x « 0.3 and Curie- Weiss behavior 
for x ~ 0.7 jl| (for further references see Q). 

These observations are counterintuitive given our con- 
ventional understanding of Na^CoO^. In NaCo02, Co is 
in the 3+ configuration and thus has six electrons. The 
cubic component of the crystal field splits the d mani- 
fold into a set of 3- fold t 2g orbitals and 2-fold e g orbitals, 
while the trigonal component will further split the t 2g or- 
bitals into dig and e' g . In this scenario, the six electrons 
of the Co will fill the t 2g orbitals (ie. a\ g + e' g ) resulting 
in a band insulator. On the contrary, in C0O2 the Co 
will be in a 4+ configuration having 5 electrons in the 
t 2g shell, and may either be a metal or a Mott insulator. 

There are two main puzzles posed by the above listed 
experimental observations. First, given that NaCo02 is 
a band insulator, it is difficult to understand why correla- 
tions are observed for the nearby composition of x = 0.7. 
The density of holes will be relatively low and therefore 
the on-site coulomb repulsion will have a minimal effect. 
Second, one would expect correlations to increase as the 
system is doped towards x = 0.3 as the hole density in 
the t 2g shell is increasing towards integer occupancy with 
an open shell. We resolve this puzzle by proposing a min- 
imal low energy Hamiltonian which captures the essential 
physics and solve it within LDA+DMFT. 
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where i, j are site indices, a, (3 are orbital indices run- 
ning over dig and e' g , t is the hopping parameter, U is the 
on-site Coulomb repulsions, and e is the on-site potential 
which mimics the Na potential. It should be emphasized 
that the orbitals to which the creation operators refer 



are low energy orbitals, effectively composed of oxygen p 
orbitals as well as Co d orbitals which compose the t 2g 
anti-bonding manifold. All the parameters of the Hamil- 
tonian are in general functions of the total density, as the 
system rehybridizes as a function of doping. 

The hopping elements of the Hamiltonian are obtained 
using LDA calculations on NaCoC>2 and C0O2 to be used 
for DMFT calculations in the Na rich region and Na poor 
region, respectively. The projected DOS for the a\ g and 
e' g states are shown in figure El and respectively. In 
both cases, the a\ g orbital has the overwhelming frac- 
tion of hole density. The doping dependence of the hop- 
ping parameters is substantial given the changes in the 
shape of the DOS and given that the total bandwidth 
increases by 0.25 eV for C0O2. Detailed cluster calcu- 
lations for LiCo02 fit the on-site U to the experimental 
photoemission yielding U = 3.5eV 0. This should be 
very similar to NaCo02, and we use a slightly reduced 
value of U = 3eV because we are considering low-energy 
orbitals. Additionally, we assume U to be independent of 
doping. A nearest neighbor Coulomb interaction should 
be added to eqnin order to assist in forming charge or- 
dered phases |5j , and in single-site DMFT this results in 
a Hartree shift. 

The on-site potential e is treated as a random binary 
variable, and the fact that it is binary will justified in 
the LDA calculations below. Treating the potential as 
random should be accurate for average quantities in the 
Na ordered phases, or even k-resolved quantities in the 
disordered region of the phase diagram. A fraction X\ 
of the lattice sites will have a on-site potential t\ while 
/jdl — x\) sites will have a have a on-site potential e 2 = 
*\-\ ■ We will assume that x ei corresponds to 1 — XNa, 
where xno, is the fraction of Na in the material. This also 
will be justified by our LDA calculations. 

We propose that the Na potential plays a critical role 
in the behavior of this material. In a previous study of 
Lij;Co02 0( a very similar material), LDA calculations 
demonstrated that in the dilute Li-vacancy region of the 
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phase diagram, the Li vacancy binds the hole which is 
usually doped into the t2 S band and hence forms a half- 
filled impurity band split off from the valence band (anal- 
ogous to Phosphorous doped Silicon, except the vacan- 
cies are highly mobile). Due to the fact that the im- 
purity band is half-filled, it may be strongly correlated 
and form a Mott insulator if the on-site interaction is 
sufficiently strong. This behavior persists until enough 
holes are present to screen the Li vacancies upon which 
the more standard picture of doping into the ti g states 
is recovered. The LDA calculations predict the impurity 
band to completely merge into the valence band once 
x w 0.98. Due to the similarities of the two materials, 
it is likely that a similar scenario occurs in Na^OoC^. 
Of course, the experiments of interest are in a region of 
much higher doping (ie. i « 0.7) where it is implausible 
that an impurity band would persist. However, it is likely 
that the Na vacancies will still act as a strong perturba- 
tion towards the Co sites. Therefore, it is possible that 
the Na vacancy potential will cause certain Co sites to 
be higher in energy and possess a higher occupation than 
the average number of holes. Increasing the occupation 
of a given site towards integer filling will favor moment 
formation assuming that the on-site coulomb repulsion is 
sufficiently strong. This phenomena has been studied in 
detail in the context of model Hubbard Hamiltonians . 
Merino et al. suggested a Hubbard model on a triangular 
lattice in the presence of an ordered potential in order to 
explain the correlated behavior observed in Na x Co02 0- 
They note that in the limit of strong coupling of the or- 
dered potential that correlations will be formed, similar 
to our work. NMR measurements strongly support this 
notion of certain Co sites being favored as there are dis- 
tinct Co sites for x = 0.7 while there are not for x = 0.3 
||. The fact that only one site is observed at x = 0.3 
supports the notion that the Na potential is well screened 
an no longer a significant perturbation for small x. 

In order to understand the effect of the Na potential 
described by the on-site potential e in eq. Q we P er ~ 
formed LDA calculations on the x = i and x — | unit 
cells which were experimentally determined by Zandber- 
gen et al0- These unit cells have 6 and 16 formula units 



for the 



and 



— | unit cells, respectively, and 
both contain equal numbers of occupied Nai and Na2 
sites. The Nai sites project directly onto the triangular 
Co lattice while the Na2 sites project onto the centroids 
of the Co triangular lattice . All LDA calculations were 
performed using the Vienna Ab-initio Simulation Code 
(VASP) [13 ■ In order to quantify the effect of the Na 
potential on the different orbitals, the first moment mj a 
of the projected DOS (ie. the on-site orbital energy) of 
the Oxygen and Cobalt orbitals arc calculated for each 
site of the supercell. m.i a = y~ L 7~^~ where the indices 
i and a arc the site and orbital indices, respectively The 
deviation of the orbital energy from the site average is 
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FIG. 1: The on-site orbital energy,m ici , relative to the average 
for various orbitals at each site in the unit cell. The energies 
for the oxygen orbitals have been divided by ten. The left 
and right panels correspond to Na3CoC>2 and Nai_Co02, re- 
spectively. 



shown in figure [T](ic m ia —fh a ). Remarkably, the results 
are roughly binary distributions as evidenced by the hor- 
izontal lines which are formed by the data points. The 
Na ordering causes each site to be perturbed, resulting in 
either a high or low energy orbital on a given site. The 
oxygen orbitals are perturbed the most followed by the 
&i g and e' g orbitals. The perturbations are substantially 
smaller for x=^ as compared to x=|, which is expected 
given that a larger amount of hole density is present at 
x=-| to screen the Na potential. A very simple rule gov- 
erns the observed behavior. In the x=| structure, all 
of the low energy ai 9 and e' g orbitals have one nearest- 
neighbor Nai present, while the high energy orbitals do 
not. There are 6 Nai sites occupied in this unit cell, each 
of which will have two Co nearest neighbors, and hence 
there are 4 high energy Co sites and 12 low energy sites. 
For the x=| structure, the ai 3 orbital follow the same 
rule. However, the e' g orbitals are perturbed in the oppo- 
site direction in this case. Although, the perturbations 
of the e' g orbitals are all less than 2meV and should be 
considered carefully. The splitting of the oxygen orbitals 
also follow a simple rule in the x=| structure. All of the 
low energy oxygens have three nearest-neighbor Na while 
the high energy oxygens only have two nearest-neighbor 
Na. For the x=| structure, the oxygen do not split into 
a binary distribution, but it can be roughly understood 
as a ternary distribution if one analyzes both nearest and 
next-nearest neighbor Na. The low energy oxygen have 
three neighboring Na, the high energy oxygens have two 
Nai neighbors, while the intermediate energy oxygens 
have two Na2 neighbors. Regardless, the a% g and e' g or- 
bitals seem to depend only on the positions of the Nai. 
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This seems very reasonable for the ai g orbitals given that 
they point directly towards the Nai sites (ie. a\ g = d z i 
in the hexagonal coordinate system). This can be very 
important for several aspects of the low-energy behavior 
in this material. For example, this effect would tend to 
drive charge ordering to occur on the Co sites which are 
located directly above/below the occupied Nai sites. Ad- 
ditionally, this could effect the stabilization the pockets 
in the Fermi surface as the two orbitals are perturbed 
differently. For a recent study on the effect of Na on the 
Fermi surface see [f]| . 

Obtaining parameters for a model Hamiltonian from 
LDA calculations is still an open problem. Therefore, 
we simply extract the qualitative effects of the Na po- 
tential. The binary splitting of the oxygen for x = § is 
roughly 350 meV while for x = ^ the maximum splitting 
is roughly 50 meV. This should be relevant to the low- 
energy Hamiltonian given that the hole density on the 
tig states has associated hole density on the oxygen via 
the rehybridization mechanism The changes on the 
Co are smaller, with the a\ g splitting being roughly 55 



and 20 meV for x = | and x = j, respectively. This is 



to be expected given that these states are near the Fermi 
energy and are not only well screened, but overscreened 
by LDA, and thus their bare values may be significantly 
larger. We take a value of ei = iOOmeV for x = | and 
ei = ^Y-meV for x = |, reflecting the fact that x — \ 
more screened. A more appropriate value for x = ^ might 
be ei = given that only a single Co is observed in NMR 
H, and we explore this as well. We conservatively as- 
sume e to be orbitally independent, which will make it 
more difficult to form correlations near the band insula- 
tor. Below we will demonstrate that these estimates yield 
predictions consistent with experimental observations. 

With these order-of-magnitude estimates, we shall now 
proceed to solve the prop osed Hamiltonian (equation 
within LDA+DMFT Ha. i ncluding the additional effect 
of binary disorder (see [l3j for references). In this case 
there will be two impurity models which represent the 
two different binary environments. We begin with a guess 
for the average hybridization function A, and then con- 
struct the bath function Go = {iuj n — Ei„ lp + /i — A) -1 
where Ei rnp is the average impurity level and /x is the 
chemical potential. We then construct the two disor- 



der bath functions 



(G 1 - ci)" 



and G 2 



(Gq 1 — £2) 1 ■ The two corresponding impurity prob- 
lems are then solved and the average Green's function is 
constructed as G = x ei G ei + (I — x ei )G e2 . The average 
self-energy is then constructed using Dyson's equation 
E(iw n ) = G^iiUn) - G _1 (io;„). Finally the DMFT 
self-consistency condition is performed: 



G (iuj n ) 



£(i^„) + (J2k iu n -H k +»-E(iu n ) 

where Hk is the homogeneous Hamiltonian (ie. eq 



(2) 



^ with e = 0) in k-space. The entire process is then 
repeated until convergence is achieved. Given that the 
local Green's function is diagonal (ie. a\ g + e' g ), the sum- 
mation over k in eqn. [3 can be replaced by a Hilbcrt 
transform (see 0). 
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FIG. 2: Spectral functions for Nao.7Co02 with random binary 
disorder for the ai s and e' g orbitals for T=290 K. 

The DMFT impurity problem is solved using Hirsch- 
Fye Quantum Monte-Carlo (see [fillip 

using the LISA 

code [fj]- All runs were performed with SR = 1 and 
500,000 Monte-Carlo sweeps. We begin by considering 
the Na rich case of x = 0.7. The spectral functions for 
the aig orbital clearly shows the correlations increasing as 
the on-site potential is increased (Figure |2J. The states 
at the Fermi energy are suppressed and spectral weight is 
transferred from low energies to higher energies as the on- 
site potential is increased. In both cases, the e' g orbitals 
are pushed beneath the Fermi energy and the pockets 
are destroyed. Another useful quantity to analyze is the 
local magnetic susceptibility (see figure|3J), which is one of 
the key experimental measurements. The susceptibility 
is relatively flat for the case with zero on-site potential, 
while a Curie- Weiss tail is clearly formed for the case 
when a on-site potential of 0.4 eV is included. In that 
latter case, 30% of the sites have 0.73 holes while 70% 
of the sites have 0.1 holes. This is qualitatively similar 
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FIG. 3: Local Magnetic susceptibility versus temperature for 
NasCoCb. 



to what is observed in NMR experiments [8j- Thus the 
ei site becomes the favored site receiving significantly 
more hole density and being driven towards integer-filling 
where a local moment forms. 

Now we proceed to the Na poor region and perform 
calculations at i = -|. The spectral function for the 
a\ g orbital shows some indication of correlations as some 
spectral weight has been transferred to higher energies 
(see figure 0J. Additionally, the e' g orbitals show a small 
occupation in this case and t he p ockets survive as pre- 
viously shown by Ishida et al |l5j . The on-site potential 
only has a small effect on the spectral function. The main 
result is that the magnetic susceptibilities only shows a 
weak temperature dependence in both cases and there 
is no sign of a Curie- Weiss tail (see figure |3J) • The on- 
site potential clearly enhances the susceptibility but is 
not capable of building strong correlations. The on-site 
potential creates two distinct sites with 70% having 0.76 
holes and 30% having 0.5 holes. Both susceptibilities 
cross the Curie tail of the x = 0.7 case in the vicinity of 
room temperature, similar to the crossover observed in 
experiment between 150-250 K [[J. 

An additional LDA+DMFT calculation for pure C0O2 
demonstrates that it is a strongly correlated metal with 
well formed Hubbard bands and a quasiparticle reso- 
nance at the Fermi energy. This guarantees that the 
susceptibility calculated for Nao.3Co02 will not develop 
a Curie- Weiss tail at low temperatures which could not 
be reached in our QMC calculations. 

In conclusion, we have proposed a low-energy Hamil- 
tonian that is capable of explaining the observation of 
correlations near a band insulator and lack of correlation 
in the Na poor region of the doped cobaltatcs. Realis- 
tic parameters for the Hamiltonian obtained from LDA 
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FIG. 4: Spectral functions at Nao.3Co02 for the a lg and e' g 
orbitals for T=290 K. 



calculations reproduce the generic behavior observed in 
the experimental magnetic susceptibility measurements. 
The Na potential is shown to be a key element in forming 
correlations near the band insulator. LDA calculations 
for the experimentally predicted Na ordcrings imply that 
the Nai sites act as a binary perturbing potential of the 
Co sites. Our model may be further improved by treat- 
ing spatial correlation of the on-site potential, improved 
downfolding techniques to obtain the parameters of the 
low-energy Hamiltonian, inclusion of polarons, and the 



inclusion of the e g states. 



The e g states are needed to 



describe the high energy behavior of this system, and 
additionally may cause further renormalization at low 
energies. Including these effects may allow for signifi- 
cantly smaller values of the on-site potential e in order 
to achieve the same qualitative effects. 
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